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1. INTRODUCTION 

In recent years, the renewable energy/green energy technologies are growing rapidly and they are 
set to reach ambitious targets as they have the ability to reduce energy crisis with environmentally friendly 
manner. Denmark has set a target of 50% of its power demand from renewables by 2030, whereas Germany 
has set a target of 35% and UK has set a target of 30% of its power demand from renewables by 2020 [1]. 
Renewable energy sources (RESs) such as wind power and solar photovoltaic (PV) energy are often spatially 
distributed and inherently intermittent and unpredictable. The power output of RESs depends on the wind 
speed, solar irradiation, and other climatic factors, which are uncertain and intermittent in nature. 
Various countries are moving towards the microgrid with RERs and other locally available distributed 
generation sources to get the benefits of providing power to remote areas with high power quality and 
reliability. Microgrid (MG) is an active distribution network, which can operate in both island and connected 
modes. In island mode, MG is disconnected from the main grid, and it supplies their loads with existing 
distributed generations (DGs). The majority of DGs used in MGs are renewable types, which are connected 
to the grid by inverters as switch mode power supplies. Microgrids with RESs and electric vehicles (EVs) can 
reduce the emissions, loss reduction, costs of energy supply, peak shaving, enhances the reliability, voltage 
profile and improves power quality [2]. 

The integration of wind, solar PV energy systems, EVs, BES systems and demand response (DR) 
provides additional energy and storage options to the microgrid. The survey on existing research works on 
the optimal strategies for the operation of MGs is proposed in [3]. Reference [4] reviews and compares 
different optimization techniques applied for achieving optimal operation of MGs from diverse perceptions, 
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in terms of objective functions, optimization techniques, constraints and computational performances. 
An operation optimization modeling for a MG considering the distributed generation, environmental factors 
and DR is proposed in [5]. Reference [6] presents an optimal power balancing model of distributed systems 
with distributed generation and EVs based on quadratic rotated conic programming. An optimal day-ahead 
scheduling and control models for MGs under uncertainty is proposed in reference [7]. A mathematical 
model to compute the optimal operation schedule that embodies DR in MG is presented in [8]. The impact of 
DR program intensity and BES system size on the operation of networked MGs is analyzed while 
considering the prevailing uncertainties is proposed in [9]. In [10], the DR program and battery are utilized to 
tackle the stochastic nature of wind, solar power generation and load demand to utilize an energy 
management system to manage hourly power of MG and optimally supply the demand. An approach to 
review the role of renewable energy as an energy producer in unit commitment scheduling problem is 
proposed in [11]. 

An optimal robust model for short-term operational scheduling of distribution network in the 
presence of uncertainties to minimize the cost of network operation is proposed in [12]. Utilization of PV 
systems for EVs charging for transportation requirements of smart cities is proposed in [13]. An economic 
dispatch problem of MG, i.e., to minimize the fuel cost of dispatchable DGs by using four different 
optimization techniques - lambda iteration, lambda logic, direct search method, and particle swarm 
optimization is proposed in [14]. In [15], energy storage units are considered as an effective option for 
handling the imbalance of wind power generation. Reference [16] proposes an approach to integrate power to 
gas technology to the multi-source MG including wind power, electric boiler, micro turbine and gas boiler. 
In [17], the energy storage systems and DR are introduced to the traditional scheduling model of wind and 
thermal power with carbon emission trading for reducing the influence of random nature of wind power on 
system operation. From the above literature review, it can be observed that there is a pressing requirement for 
an optimal operation of MG considering the RESs, BES systems, EVs and DR [18-20]. 

The main contribution of this paper is to optimize the total operating cost of microgrid (MG) by 
considering the power exchange between main power grid and MG, renewable energy sources, battery 
energy storage system, electric vehicles and demand response. In this paper, an optimal scheduling 
methodology is presented for the MGs in large distribution systems. Here, the energy storage batteries and 
EVs are used to improve the operational flexibility of the system and to increase the penetration of RESs. 
The demand response that is considered in this paper will work as an operating tool for the improvement of 
reliability and efficiency of MG. The aim of demand response is to encourage the customers to become active 
partners of MG energy management system. 


2. MODELING OF MICROGRID COMPONENTS 

Modeling of MG components such as WEGs, solar PV units, load demands, BES systems and EVs 
are presented in this paper. The detailed modeling of wind energy generators, solar PV units and load 
demands and their uncertainty handling approaches has been presented in references [21-23]. 


2.1. Modeling of Battery Energy Storage (BES) System 

Generally, the battery is considered as constant voltage source. BES system is used to smooth out 
the energy supply by absorbing energy when the power generated by renewables is more than the scheduled 
power, and dispensing energy when the renewables do not fulfill the load demand. Charge/discharge power 
limits of battery are expressed as [24-25], 


ch,t chmax 
C2P ah, (1) 
disch,t disch,max 
gis paseht = pat (2) 
where Pf?’ and PZis’™* are the maximum charging and discharging rates of battery [26]. The state of 


charge (SOC) of battery plays a major role in controlling the charging/discharging and in optimal scheduling 
of battery [27]. The constraints on SOC of storage battery are expressed as, 


SOC eS S00; = SOC (3) 
t CBEsk 
sock = (7) ” 
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where SOC?" and SOCj"™ are minimum and maximum BES capacity limits, respectively. The C£es, can 
be expressed as [28], 


ht disch,t 
Chase = Cine + (F re nen) — ee Naiscn) (5) 
where C;,,;, is the initial capacity of battery storage system. 


2.2. Modeling of Electric Vehicles (EVs) 

EVs will play two roles when they plug into a MG. These EVs will appear as loads during the 
charging, and they also used as storage units to supply energy to the remaining load demands, if the batteries 
have the sufficient spare energy. The battery of EV acts as a source of stored energy to provide a number of 
ancillary services to supply primary and secondary control and voltage control. The battery systems of EV 
are the chemical storage devices, and they are exponential functions over time. The instantaneous charging 
status of the battery system of EV is expressed as [29], 


PE, = PHHE*(1 — ett/tmax) + Po, (6) 


where Pg, is the current status of battery system, P72** is the maximum power capacity of EV, tym, is the 
maximum charging time. Similarly, the discharging status of the battery system of EV is 
expressed as [30], 


Pay = Pyy.e{t/tmaz) (7) 


where a@ is a constant parameter and it is calculated by considering a fully empty battery system of EV that 
absorbs 97% of maximum power capacity, and its value is approximately equal to (tmg,/3)- 


3. PROBLEM FORMULATION 

The proposed structure of MG that has been considered in this paper is depicted in Figure 1. In this 
paper, the power exchange between main power grid and MG is considered. The components that are 
considered in the proposed MG structure are: distributed generators (i.e., WEGs and solar PV units), 
BES systems, EVs and demand response. The EV aggregator collects the information of 
charging/discharging strategies of EV batteries in electric vehicles [31]. The EV aggregator participates in 
the market clearing process to reduce the total operating cost (TOC) of the system by determining the mode 
of operation of EV and to provide an operational flexibility. 
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Figure |. Proposed Structure of Microgrid (MG) 
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The objective considered in this paper is to minimize the total operating cost (TOC) of the system, 
which is the sum of cost due to power exchange between main power grid and MG, cost due to power output 
from WEGs, solar PV units, BES systems, EVs and demand response [32]. Mathematically, it can be 
expressed as, 

Minimize TOC, 


TOC = YhealCéria Paria) + Diet Civ (Piva) + Djer C$ (Phv,j) + Lees C5(P5. x) + DY Cv (Piv.) + Cor Pon) (8) 
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The first component in (8) represents the cost due to power exchanged between main grid and MG, 
and it can be expressed as [33], 


t t jy=ct t 
Coria (Péria) = Coria PGria (9) 
where Pé,,q is the power output from main power grid at hour t. Pé,;q is positive if the main grid transmits 


power to the MG, and it is negative if the main grid absorbs power from the MG. Second term in (8) 
represents the direct cost associated with wind power and it is expressed as [34], 


Ci (Pi) = Cb Pe (10) 


Third term in (8) represents the direct cost for the scheduled solar PV power, and it is expressed as, 
CE (Ply) = C8. Phy 4 (11) 


Fourth component in (8) represents the operating cost of BES system, which includes the 
installation and degradation costs, and it can be expressed as, 


Ce St eee (12) 


Fifth component in (8) represents the cost associated with the charge/discharge power of EV, and it 
can be expressed as, 


Civ (P41) = Chy Phy (13) 


Last component in (8) represents the cost due to incentive based demand response (DR), and it can 
be expressed as, 


Chr (Por) = Cor: Por (14) 


where P§, is the amount of planned power/load reduction, Ch, is the incentive payment for load reduction. 
The above objective function is solved subjected to the following equality and inequality 
constraints. 


3.1. Power Balance Constraints 
This equality constraint of MG system maintains the balance between power generation/supply and 
demand. This constraint is expressed as [35], 


N N N N N 
Porid + pas Pw i + Der Pov j + iar Pik + biel Pvt + Dea EDR = Pr + Pies (15) 


where t = 1,2,..., T. P} is the total system load at hour t. 
3.2. Power Generation Capacity Constraints 


For stable operation of MG, the active power output of each MG component is limited by lower and 
upper bounds, and they are expressed as [36], 


ata at epics ve (16) 
Pw i < Pwr (17) 
Poy, S Poy (18) 
Pha & Phrmax (19) 


where Py, ; is the power generated by ith WEG at hour t, Py'{* is its upper limit, Psy, j 18 the power 
generated by the jth solar PV unit at hour t, Pay; is its upper limit, and P5p.max is the maximum allowed 


load reduction at hour t. 
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3.3. Constraints on BES System 
The state of charge (SOC) balance equation of a BES system is expressed as, 


Soct = Soct1 + I‘At (20) 


where I* is the charge (+)/discharge (-) current (in kA). The charge/discharge current limits are expressed as 
[26], 


pmin < qt < pmax (21) 


3.4. Constraints on Transmission Line Capacity 
As mentioned earlier, the capacity of transmission line that connects between the main power grid 
and MG can be expressed as, 


Sprite Peg pee (22) 


where P/"** is the maximum capacity of transmission line. 


4. RESULTS AND DISCUSSION 

The Microgrid (MG) considered in this paper has the WEGs, solar PV units, battery banks, electric 
vehicles (EVs), demand response and the power exchange between MG and main power grid. In this system, 
four wind mills with the capacity of 750kW each, four solar PV plants with the capacity of 250kW each, 
two battery units with the capacity of 250kW each and two electric vehicle aggregators with the capacity of 
200kW each are considered. Hourly load demand, wind speed and solar irradiation data for 24 hour period 
has been presented in Table |. The optimization problem formulated in this paper has been implemented by 
using Generalized Algebraic Modeling Systems (GAMS) software and solved by SNOPT solver [36]. 
All the developed programs were run on a personal computer with Intel Core i7, 3.0GHz and 8GB RAM. 
In this paper, V,;, Veg and v, are considered as 3m/s, 25m/s and 10m/s, respectively. 


Table 1. Load Demand, Wind Speed and Solar Irradiation Data for 24 Hour Period 
8 Load Wind speed Solar irradiation Hour Load Wind speed Solar irradiation 








(kW) (m/s) (W/m?) (kW) (m/s) (W/m?) 

1 2200 5.7 0 13 5300 5.9 833 
2 2400 6.5 0 14 5200 4.9 850 
3 2800 Bs) 0 15 5000 3.5 680 
4 3200 6.9 0 16 4600 3.4 595 
5 3400 8.6 93.5 17 4000 2.8 255 
6 3800 10.5 212.5 18 3700 3.1 212.5 
7 4000 13.6 255 19 3400 2.3 153 
8 4200 10.4 467.5 20 3200 2.9 68 
9 4600 9.1 637.5 21 3000 3.5 42.5 
10 5000 9.3 680 22 2800 3.8 0 
11 5200 V7 816 23 2600 3.8 0 
12 5400 7.0 850 24 2400 4.8 0 





In this paper, four different case studies are performed, and they are described next: 
a) Case 1: Considering the power output only from main power grid. 
b) Case 2: Considering the power output from main power grid and renewable energy sources (RESs). 
C) Case 3: Considering the power output from main power grid, RESs, BES systems and electric 
vehicles (EVs). 
D) Case 4: Considering the power output from main power grid, RESs, BES systems, EVs and demand 
response. 


4.1. Case 1 

In this case, the objective function consists only first term of (8), (i.e., cost due to power output from 
main power grid. As mentioned earlier, this objective function is solved by using the SNOPT solver of 
GAMS software, and the obtained results are reported in Table 2. The total operating cost (TOC) obtained in 
this case is 590158.5$/day, and this cost is due to the power output from main power grid. 
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4.2. Case 2 

In this case, the objective function consists only first three of (8), (i.e., cost due to power exchange 
from main power grid and MG, and renewable energy sources (1.e., wind farms and solar PV plants)). 
The forecasted and actual power outputs of WEGs and solar PV units are depicted in Figure 2. Table 2 also 
presents the optimum objective function values for this case. The optimum cost incurred in this case for the 
24 hours scheduling period is 515792.1$/day, which is the sum of cost of power exchange between main 
power grid and MG (ie., 116604.3$/day), cost of wind power generation (301567.2$/day) and the cost of 
solar PV power generation (97620.6$/day). The total operating cost (TOC) obtained in this case is 12.6% less 
than the cost obtained in Case 1. 





Forecasted and actual wind and solar PV powers for Case 2 
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—*— Forecasted wind power 
—S— Actual wind power 

—*— Forecasted solar PV power | 4 
—S— Actual solar PV power 


2500 - 











2000 


1500 


1000 


Forecasted and actual powers (in kW 
a 
° 
[o) 











L L oi» See @-¢ 
10 25 
Time (in Hours) 


Figure 2. Forecasted and actual wind and solar PV powers for Case 2 


Table 2. Optimum Objective Function Values for Cases 1, 2, 3 And 4. 





Objective Function Values Case 1 Case 2 Case 3 Case 4 

Cost of power from main power grid (in $/day) 590158.5  116604.3 71512.6 45762.5 
Cost of wind power generation (in $/day) — 301567.2  295391.3 284691.8 
Cost of solar PV power generation (in $/day) — 97620.6 91359.7 87368.4 
Cost of battery storage and EVs (in $/day) — — 42128.2 39850.5 
Cost of demand response (in $/day) = — — 36484.2 
Total operating cost (in $/day) 590158.5 515792.1 500391.8 494157.4 

Cost reduction (in %) — 12.6 15.2 16.3 





4.3. Case 3 

In this case, the objective function consists all the terms of (8) except last term, (i.e., cost due to 
power output from main power grid and MG, renewable energy sources (i.e., wind farms and solar PV 
plants) and battery energy storage (BES) system). In this case, to extend the lifetime of BES system, the 
soci” and SOC} of the BES system are set to be 20% and 80%, respectively. This BES system stores 
the energy during the off-peak hours by charging the battery until it reaches SOC{"**; whereas during the 
peak hours, it injects the stored energy to loads until it reaches SOC7™”. 

The TOC obtained in this case is 500391.8$/day, which is the sum of cost due to power exchange 
between main power grid and MG (ie., 71512.6$/day), cost due to wind power generation (i.e., 
295391.3$/hr), cost due to solar PV power generation (i.e., 91359.7$/day) and the cost due to BES system 
and EVs (i.e., 42128.2$/day). From the simulation results, it can be observed that the TOC obtained in this 
case is 15.2% less than Case 1, and 3% less than Case 2. 


4.4. Case 4 

In this case, the objective function (i.e., TOC of the system) consists all the terms of (8), (i.e., cost 
due to power output between main power grid and MG, renewable energy sources (i.e., wind farms and solar 
PV plants), BES system, EVs and demand response). The TOC obtained in this case is 494157.4$/day, which 
is 16.3% less than Case 1, 4.2% less than Case 2, and 1.3% less than Case 3. 
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From the simulation results of four case studies, it can be observed that by operating the MG with 
RERs, BES system, EV and demand response simultaneously can reduce the total operating cost of the 
system, which leads to the optimal utilization of hybrid energy sources in the microgrid. 


5. CONCLUSION 

This paper proposes an optimal scheduling approach of Microgrid by optimally utilizing the 
renewable sources of energy such as wind and solar photovoltaic (PV) powers. The optimal operation of 
Microgrid is implemented by including the stochastic nature of wind energy generators, solar PV systems, 
battery storages, electric vehicles and optimal dispatch of schedulable loads, i.e., demand response. In order 
to reduce the intermittent nature of wind, solar PV powers, electric vehicles (EVs), battery energy storage 
(BES) systems and demand response are introduced in this paper. Minimization of total operating cost (TOC) 
(i.e., cost of power exchanged between main power grid and MG, wind, solar PV powers, BES systems, 
EVs and demand response) of the system is considered as an objective function. In this paper, the probability 
distribution functions are used to handle the uncertainties related to wind power, solar PV power and load 
demand. The proposed optimization problem has been solved using the GAMS software. The simulation 
results of this paper show that the total operating cost of the system has been reduced significantly by 
optimally utilizing various microgrid components simultaneously. 
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